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, A lattice boson model is used to study ordering phenomena in regular 2D array of supercon- 

' ductive mesoscopic granules, Josephson junctions or pores filled with a superfluid helium. Phase 

diagram of the system, when quantum fluctuations of both the phase and local superfluid density 
are essential, is analyzed both analytically and by quantum Monte Carlo technique. For the system 
of strongly interacting bosons it is found that as the boson density no is increased the boundary of 
ordered superconducting state shifts to lower temperatures and at no > 8 approaches its limiting 
position corresponding to negligible relative fluctuations of moduli of the order parameter (as in an 
£f~) ' array of "macroscopic" granules). In the region of weak quantum fluctuations of phases mesoscopic 

phenomena manifest themselves up to no ~ 10. The mean field theory and functional integral 1/no 
- expansion results are shown to agree with that of quantum Monte Carlo calculations of the boson 
Hubbard model and its quasiclassical limit, the quantum XY model. 
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I. INTRODUCTION 



The study of mesoscopic systems has resulted thus far in many new interesting fundamental concepts . Progress 
has been especially rapid due to the development of nanotechnology methods which has opened up new avenues for 
sophisticated experiments. In this connection the study of arrays of ultarsmall granules, microclusters or Josephson 
junctions is of particular concern (see e.g. 
1 I Granular superconductors, Josephson arrays, superfluid helium in a porous media Q 3X6 . 3S cL rule, described in 
terms of different modifications of quantum XY model (see below) , but this description is correct only if relative fluc- 
tuations of the local superfluid density are not essential M . It takes place in sufficiently large granules at temperatures 
far below that of onset of superconductivity in each individual granule. To study the role of quantum fluctuations of 
moduli other more adequate models should be used. 

A convenient starting point for the description of the N x N system of interacting bosons (Cooper pairs in granules, 
r ^ ' He atoms in pores etc.) is the Bose - Hubbard Hamiltonian: 

cn ' t U 2 

Tj" \ flh = 7 V {ZoJiCLi - a^idj - a ] jdi} + — {ahai - n } (1) 

■ <hj> 1 

00 

where a] (dj) is a boson creation (annihilation) operator at a site i = 1, N 2 ; t is the strength of the hopping between 
nearest neighbor sites < i,j > and U > is an on - site repulsive interaction. 

The system (|l|) has a rich phase diagram |||]], containing a Mott insulating phase (at zero temperature), the 
superfluid and normal (metal) phases. At a commensurate density no =< a\a,i > (number of bosons is an integer 
multiple of the number of sites) and T = the boson Hubbard model lies in the same universality class (see [|| - Q ) 
as the quantum XY model with Hamiltonian: 

H xy = J J2 i 1 - cos to - <Pi)} - t E i d / d ^} 2 (2) 



o 
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where tp^ S [0, 2ir) are phases of the order parameter. Obviously, at T 7^ the requirements the density to be 
commensurate uq — k is too stringent. In the latter case, the behaviour of the system (|l|) will depends continuously 
; 1 1 on no j and the critical properties will be the same in some band uq = k ± Sno, the width 2<5no should decrease as the 
temperature is lowered. 

The properties of the system (^), that have at finite temperatures the superfluid and metallic phases are described 
by two dimensionless parameters: the temperature in units of Josephson coupling constant T — ht,T/ J and quantum 
parameter q = y/U/J which is responsible for the strength of zero - point fluctuations of phase. Corresponding 
parameters of the Hubbard model are T — k^T j (trio) an d 1 = \jU /(trio). 

The general purpose of this communication is comparing of phase diagrams of models ([!]) and (^|) to estimate 
the importance of the mesoscopic phenomena in regular 2D systems. Of prime interest to us is the case of finite 
temperatures, apart of intriguing quantum phase transitions which take place at T = (see e.g JTT[ ] and references 
therein). Sections || and III present the mean field and functional integral 1/no - expansion approaches. From ab 
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initio quantum Monte Carlo calculations of different characteristic quantities (Section IV) we determine the phase 
diagram T£(q;no) of the boson Hubbard model at different densities no and compare it (Section [v]) with the phase 
diagram of 2D quantum XY model. 

II. MEAN FIELD APPROXIMATION. 

A qualitative estimation of the phase diagram of the boson model (0) can be obtained in a simple mean field 
approximation (see e.g. |0|, Jl^,[l4| and references therein). The boundary T c (q\ no) of ordered state can be found in 
MFA from the equation: 

2 E {n + n + 1} L-i 2 [n-n? /[2T\ _ e - q -[n+i-^/[2T] j /{2n + i _ 2 n} 

1 _ n=-n ' ,„s 

~ _ 53 v>) 

Z n YL e-? 2 {"-"} 2 /{2T} 
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where n = fi/U — z/(2q 2 no), z is the number of nearest neigbours (z = 4 for 2D square lattice). The Equation ([?]) on 
the boundary of ordered state differs in some details of that of the Reference |l3| . Both of them became equivalent 
in the limit of large densities n , with the Equation (^) more accurate at n ~ 1. 

The condition on the chemical potential is that the mean number of particles be equal to no: 

oo 

ncxp(-q 2 {n-r 1 } 2 /{2T}) = Q (4) 

n— — no 

It should be pointed out, that in the limit no — ► oo Equation (|j) gives u = and the boundary of ordered state 
transforms to that of quantum XY model |Q . 

The lines T£(q;no) of the boson Hubbard model obtained from Equations ([|) and (Q) are shown in Figure 1 for 
different densities no. Calculation shows that the lines T£(q;no) reaches their limiting position, the phase boundary 
T£ y (q) of quantum XY model at no > 25. 

III. 1/N - EXPANSION. 

To improve the qualitative mean field estimation of the difference of phase diagrams of models (|f|) and (||) let us 
represent the partition function of the Hamiltonian ([I]) in a path integral representation as a trace over a complex c- 
number Bose field <f> []lo|| : 

Z h = tr {e- s } = J £>($, (5) 

( 

$(*,*•) = j + | e i^-^i 2 + |E[i^i 2 -H ; 

[ * <i,j> i 

$ i (0) = * i C9), $*(0) = $*(/?) 
The substitution $i = + Sriie 1Vi (at integer no) in Equation (g) gives: 
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Srii = ni - n = Sriiir), w = yj»(r). 
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From Equation (g) one can see, that increasing of the mean number of particles at each granule provided that 
J = tno and U are constant, enables one to leave an action in terms of the phase degrees of freedom alone @,[n| to 
lead to the action of quantum XY model. 

Being interested in the difference between phase boundaries of models (0) and (Q) at sufficiently high but finite no, 
let us expand the superfluid density in powers of I/uq up to a second order. Defining the superfluid density from the 
response of a system to the shift of phases at the boundary [^5[ , from Equation (Q) we have: 



where 7 is the helicity modulus of quantum XY model [|16 17 1. The first order corrections are equal to zero due to 
the invariance of the action of the XY model against the " time" inversion. 

Rather a complex expression for can be represented as some equilibrium value of quantum XY model [O], 
which can be easily estimated via. quantum Monte Carlo technique or different self - consistent apprximations |p.q| . 

Given a value of the coefficient of the expansion (Q) as a function of control parameters = T^ 2 \q, T), one 
can construct an upper estimation of the phase boundary T^(q; uq) of the boson Hubbard model ([!]): 

T h ( q ,no) < T xy((? ) jl + ^—^ 1 - T xy (,) |1+ HT ^ {q) I (8) 

where T xy (g) is the line of topological phase transitions of quantum XY model (||). The estimation (||) can be easily 
obtained from the assumption that lines of phase transitions of both models are defined by the universal relation |l6| : 

7 (<?, T x c y ) = 2T x c y /7r, v s {q, T£) = 2T£/ir 

Results of above mentioned estimations are given in Figure 1. It turns out that in the region 0.7 < q < 1.5 the line of 
phase transitions of the Hubbard model approaches (to within 5%) its limit at no — 8, whereas some greater densities 
no > 16 are required in the strong quantum region q > 1.7 because of the rapid increase of the coefficient T^ 2 \q). 
Direct Monte Carlo calculation of the phase diagram T£(q) (see below), being in agreement with the predictions of 
l/no - expansion at 0.7 < q < 1.5, shows that in the case of strongly interacting system (at q > 1.7) theoretical 
estimations markedly overestimate the maximum boson density no at which mesoscopic effects are still essential. 

To conclude of this chapter one additional feature should be recognized. It is easy to see, that all estimations 
presented in this Section and started from the partition function (|5|) have been carried out in the grand canonical 
ensemble with zero chemical potential. This has enabled us to make all calculations analytically, disregarding the 
restrictions on the total number of particles in the system when integrating over the fluctuations of moduli of the 
order parameter. In order to justify the possibility of comparing the MC results with theoretical estimations it is need 
to show that, being calculated within the same approach, the discrepancy 5n between the mean number of bosons 
per granule and no is small. 

From @ for Sn one can write 



Sn= — A« + ... (9) 



no 



As calculation shows, the value of is less than 0.1 in the region q > 0.7, T < 1. This observation justifies our 
use of great canonical ensemble with zero chemical potential in estimating the coefficient of series (R). 



IV. QUANTUM MONTE CARLO SIMULATIONS. 

In studies of properties of the boson Hubbard model (||) at different values of control parameters {q, T} we have 
applied the " checkerboard" breakup [[[9) . In this method the classical degrees of freedom to be sampled are imaginary 
timedependent boson occupation number field {n? }, i = 1 . . . ./V 2 , p = . . . 4P. The algorithm of Monte Carlo (MC) 



calculations of the quantum XY model was described in 17 2C||. We have performed extensive tests to verify that 



our results converge in the limit of P — > 00. Results presented below have been obtained by the averaging over 3 — 5 
initial configurations formed by 4P - times (P - times for the XY model) multiplication of the random configuration 
of bosons (phases for the XY model) at the lattice N x N. 

The main attention has been given to the calculation of the superfluid density v s . This quantity have been 
determined from the winding number P|l9|] 
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v s = 0.hT{Wl + W*)^ (10) 

4P N 4P N 

w * = E E ^ y = E E (-i) iB+ x, 

p=0i y — 1 p— i x — 1 

where n? means the number of bosons at a site z (with coordinates {i x ,iy}) of a level p of 3D classical system. We 
have also used the current autocorrelation function pl[ 

1 ^-^rpE^^^^y (id 



n iV 2 \ x /h nlN 2 TP 



= -| E { flt ^ a » - a W,} - ^V) = e-MPjMe-r&l 1 



The substitution a, — > y/no~e ltpi transforms Equation ( |11| ) to the well - known expression for the helicity modulus 7 of 
the quantum XY model ]l7j . 

As have been pointed out by Scalapino et. al. the temperature derivative of the superfluid density gives an 
additional information about the type of phase transition at some temperature T c (q): in framework of the Kosterlitz 
- Thouless picture the value of d{[3v s )/d(3 scales to a Dirac delta function S(T — T c ). On the finite lattices d{(3v s )/d(3 
shows a response which increases with lattice size N, the position of the maximum of the derivative being independent 
on N. 

To find the derivative of the superfluid density d(Piy s )/dp we have estimated the difference in internal energies of 
systems which differ by a phase twist Sip in the boundary condition along one lattice direction. 

— ^ " S+(3 W (12) 

One can show, that for the Cooper pairs of charge 2e this phase twist can be realized in the "flux quantization" 
scheme and is equivalent to threading a flux through the center of a thorus on which the system lies [p2| . 
We have also calculated he fluctuation of bosons at lattice sites 



^ = i(EEK-o} 2 ) as) 

\ p=0 i I h 



V. RESULTS AND DISCUSSION. 



Shown at Figure 2 are dependencies of the superfluid fraction v a [T) of the Hubbard model at q = 0.2 (in the 
classical region of XY model (§), 

Figure 2a) and q = 2.0 (see Figure 2b). For reference, the helicity modulus 7 of 
quantum XY model as function of temperature T is also plotted. Analysis of data obtained at different sizes N and 
densities no of the system reveals that for the system of strongly interacting bosons (at q = 2.0) the MC results are 
in qualitative agreement with the theoretical estimations of Sections || and [II . Really, from Figures 1 and 2 one can 
see that as the density of bosons no is increased, the boundary of ordered superconducting state of the system (|l|) 
approaches that of quantum XY model with critical temperatures of the Hubbard model being greater than T^ y 
of quantum XY model. The line of transitions T c (q) can be estimated from the universal relation v s {T c ) — 2T c /tt. 
Thus defined, the temperature of metal - superconductor transition agree fairly well with the position of the peak of 
the temperature derivative of superfluid density ([l2|). Our calculation shows, that the position of the maximum of the 
derivative does not depend (to within statistical errors) on the system size and, as one can see from Figure 3, lowers 
as no is increased. 

From Figure 2a one can show, that the transition temperature of the weakly interacting bosons (at q = 0.2) 
appeared to be less than that of qua ntu m XY model. This tendency persists with increasing the system size. The 
theoretical approach used in Section III can not help to elucidate the reason of this phenomena because at q < 0.4 
relative fluctuations of moduli of the order parameter are only weakly damped by the interaction, corrections 
and A^ 1 ) are large and theoretical estimations work badly. 
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Let us consider the behaviour of v s {q) and 7(g) as functions of quantum parameter q at T = 0.5 (see Figure 4). 
Defining the point q c of phase transition from the universal relation v s (q) = 2T/n we see that the boundary of ordered 
superconducting state of the model (0) lies to the right of the boundary of the model (||). This conclusion is verified by 
the results of calculations of the derivative d((3v s )/dj3 presented at the insert of Figure 4. The position (?|n =3 ~ 2.3 
of a peak of the derivative is in a sufficiently good agreement with the critical point q c \ no =3 ~ 2.4 determined from 
the universal relation. 

The dependence of the relative fluctuation of particle number at sites of the system vs. quantum parameter q is 
shown at Figure 5. Particularly, Figure 5 can serve as an illustration of the role of interaction in the transition to 
the quasiclassical limit from the boson Hubbard to quantum XY model. Really, at finite densities no the spectrum 
of the operator hi — no can be considered as unbounded only if relative fluctuations of particle number are small 
Sn^/riQ -C 1. Then, as it is usually done on examination Josephson or granular systems in terms of the model (|^), 
the particle number operator hi — no can be chosen as a conjugate one to the "phase" operator (pi. hi — no = ld/dtfi 
Increasing of the interaction (quantum parameter q) leads to the suppressing of the relative fluctuations of the order 
parameter module as can be seen from Figure 5. It should be noted, that at high q the fluctuations of particles number 
are greater than those of quantum XY model and approach their with increase in density. The results presented in 
the insert of the Figure are relative fluctuations Sn^/n^ as functions of quantum parameter q at T — 0.5 at different 
densities uq. The increase in no is seen to be of great importance in suppressing the relative fluctuations. 

A great attention has been given thus far to the possibility of reentrance phenomena, when the global superconduct- 
ing state in some region of q is absent not only at high, but also at sufficiently low temperatures. In the framework 
of XY model the possibility of this phenomenon taking place has been connected with the domain of phases p2| , 
dissipation or mutual capacitances effects [ p3|]2~4| . From results presented one can see that taking into account the 
fluctuation of moduli of the order parameter does not lead to the reentrance phenomena at least in the region explored. 

To conclude, we have used the boson Hubbard model to analyze the effect of quantum fluctuations of phases and 
moduli of the order parameter on the onset of superconductivity in 2D mesoscopic Josephson system. Both mean 
field approximation and I /no expansion leads to the conclusion that the line T c (q) of superconductor - metal phase 
transitions lies above that of quantum XY model, the latter being quasiclassical limit (as no — > 00, U ^ 0) of the 
Hubbard one. Our MC simulations show that in the region q < 1 of small quantum fluctuations of phases it needs 
an average of 10 bosons per site to suppress relative fluctuations of the local superfluid density. As interaction is 
increased, the quasiclassical limit is approached at lower densities (no ~ 8 at q ~ 2). No reentrance or discontinuity 
phenomena have been found. 

This work was supported by Russian Foundation of Basic Research, Programs " Solid state Nanostructures" and 
Superconductivity. 
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Figure 1. 

The phase diagram of 2D Hubbard ([]]) and quantum XY (||) models. S - supercondudting, N - normal state. The 
mean - field results: 1: uq = 1, 2: no = 2, 3: hq = 6, 4: quantum XY model (no = oo); The 1/uq expansion (||) 
results: 5: uq = 6, 6: no = 14. Symbols present points of phase transitions have been found by MC method. 

Figure 2a. 

The superfluid density v s (helicity modulus 7) vs. temperature T at q — 0.2. The dependence 2T/ir (see in text) is 
given by a dashed line. Data are connected to guide the eyes. If not presented, error bars are within the size of the 
data point. 

Figure 2b. 

The superfluid density v s (helicity modulus 7) vs. temperature T at q — 2.0. 
Figure 3. 

T(3{E(-k /2) — E(0)}/{noN 2 } as a function of temperature T. The solid lines are spline fits to the data. 
Figure 4. 

The superfluid density v s (helicity modulus 7) vs. quantum parameter q at T = 0.5. The line 1/tt is shown with the 
help of a dashed line. 

Insert: the temperature derivative of the superfluid fraction (fl2"|). 
Figure 5. 

Fluctuation 5n^ of the number of bosons at sites of the system as a function of quantum parameter q at T = 0.5. 
Insert: relative fluctuations Sn^/riQ at different densities no- 
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